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1 Nuclear Magnetic Resonance QIP 



1.1 Introduction 

Quantum information science is an interdisciplinary area of study combining 
computer science, physics, mathematics and engineering, the main aim of 
which is to perform exponentially faster computation using systems which 
are governed by the laws of quantum mechanics. The task of performing 
quantum algorithms or simulations is referred to as a quantum information 
processing (QIP). More fundamentally, quantum information science pro- 
vides a new lens through which to view quantum physics. An interesting 
consequence of this new perspective is that any system which possesses an 
accurate description in terms of quantum mechanics may be simulated by any 
other quantum system of similar size [1]. Such simulations provide a link be- 
tween chemistry and quantum information science, through the discipline of 
quantum chemistry. 

In order to accurately predict the results of an atomic or molecular in- 
teraction, it is sometimes necessary to formulate a Hamiltonian model for 
the molecules in question. The Hamiltonian assigns energy values to certain 
states of the physical system, called basis states. This model is then incorpo- 
rated into an appropriate equation of motion; either the Schrodinger equation 
(for the non-relativistic case) or the Dirac equation (for the relativistic case). 
The number of basis states grows exponentially with the physical size of the 
system, making the solution of the equation of motion a difficult task for a 
classical computer. 

These interactions can be efficiently simulated on a quantum computer, 
however, since a quantum operation is effectively performed in parallel across 
all basis states taken into a superposition. Quantum simulations of the type 
described in this chapter are thought to comprise a set of attainable mile- 
stones for quantum information processing, given that there exist relatively 
simple quantum systems whose dynamics are not easily simulated using clas- 
sical computers (for example, the 10-body Schrodinger equation for electrons 
in a water molecule [2j. 

The remainder of this chapter is organized as follows: We first detail a 
series of quantum algorithms which are useful in simulation of chemical phe- 
nomena, then describe a class of QIP implementations using nuclear mag- 
netic resonance (NMR) and electron spin resonance (ESR). In conclusion, 
we discuss some recent experiments and progress toward scalable spin-based 
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implementations. 



1.2 Quantum Algorithms for Chemistry 

The ostensible goal of a chemical simulation is to extract a small amount of 
data about a given process; a rate constant, ground state energy, or other 
quantity of interest. However, it is often necessary to manipulate a large 
data structure, such as a full molecular electronic configuration, in order to 
obtain the desired output. This is the exact class of mathematical prob- 
lem for which a quantum computer is thought to be superior [U E] , often 
qualified with the notion that using a quantum system to simulate another 
quantum system provides a physically elegant intuition. In this section, we 
describe how quantum information processing can provide the desired data, 
without an exponential growth in the required computational resources. The 
algorithms used include straightforward digital simulation of Hamiltonian 
mechanics using the Trotter expansion, as well as implementations of the 
adiabatic algorithm to simulate ground state properties of a large class of 
physical Hamiltonians. These comprise the vast majority of proposed simu- 
lation techniques; though they are far from a complete list of QIP architec- 
tures. 

1.2.1 Digital Quantum Simulation 

In order to simulate continuous degrees of freedom, such as position and mo- 
mentum, it is often necessary to discretize these degrees of freedom onto a 
finite space, to ensure that the amount of memory required to store their val- 
ues is bounded, while also ensuring that the resolution remains high enough 
to mimic the dynamics of the continuous- variable system under considera- 
tion. This method was pioneered by Zalka [1| and Wiesner [3], to simulate 
Hamiltonians of the form 



To accomplish this, position is discretized as detailed above, and encoded 
into a quantum register. In this way, an arbitrary superposition of 2™ x- 
eigenstates can be stored in m qubits. 

To simulate the evolution under this Hamiltonian, it is useful to perform 
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a Trotter decomposition to first order [B]: 
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This approximation, valid for small t, expresses the evolution under the 
Hamiltonian to be simulated in terms of operators which are diagonalized 
in known bases (momentum and position, respectively). Since V depends 
only on a;, it is diagonal in the basis selected above. Assuming that diago- 
nal operators can be implemented quickly on a quantum computer, all that 
remains is to find an efficient quantum circuit that will transform operators 
which are diagonal in the x-basis into operators which are diagonal in the p- 
basis. The quantum Fourier transform fulfils this requirement, and has been 
widely studied [3 [S] and implemented [HI [IDl HI]- The effective evolution 
operator, for a single particle in one dimension, is then: 



QFT^ 
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Q FT exp [-iV{x)t], 



(3) 



where Dp2 is the diagonal operator whose entries are the eigenvalues of p"^. 

Using n m-qubit registers, one can simulate a system of n interacting 
degrees of freedom. This simulation involves the class of Hamiltonians 



H 



2m 



(4) 



Here, the index j lists the degrees of freedom in the Hamiltonian, which may 
correspond to distinct particles or dimensions. The Trotter expansion of the 
evolution under this Hamiltonian requires one quantum Fourier transform 
per degree of freedom. 

In order to benefit from the efficiency offered by quantum simulation, it 
is also necessary to ensure that measurement can be achieved in polynomial 
time. Kassal et al. [12j demonstrated that reaction probabilities, state-to- 
state transition probabilities and rate constants can be extracted from a 
quantum simulation of a chemical process in polynomial time. To accomplish 
these measurements, the classical algorithm for generating a transition state 
dividing surface is invoked, subdividing the simulation space into regions 
corresponding to products and reactants. Measurement of the single bit 
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corresponding to the presence of the wavefunction in the reactant or product 
regions can be used in conjunction with phase estimation [13], so that the 
precision of the transition probabihty scales as N~'^, where N is the number 
of single-bit measurements. 

Using a limited capacity for quantum control, it is possible to simulate 
quantum dynamics, with immediate applications to small numbers of react- 
ing atoms. Note that no Born-Oppenheimer approximation has been used. 
Indeed, the quantum algorithm is more efficient without this approximation, 
because the Born-Oppenheimer approximation requires the calculation of po- 
tential energy surfaces at many points throughout a given simulation. The 
realization of chemical simulations by the algorithm above, while scalable, 
requires hundreds of qubits in order to simulate a few particles with sufficient 
spatial precision, and ~ 10^^ elementary quantum gates to obtain sufficient 
precision in time. 

1.2.2 Adiabatic Quantum Simulation 

Given that the number of quantum gates used in the simulation detailed 
above is prohibitive for implementation in the near future, it is preferable to 
search for alternatives to gate-based quantum computation for the purpose 
of simulating certain chemical processes. One such alternative is adiabatic 
quantum computation [T3]. This method of quantum information processing 
is naturally suited to problems which relate to ground states of complicated 
Hamiltonians, for example, the problem of protein folding |15j . 

The central mechanism of adiabatic quantum computation is the transfor- 
mation of a prepared state to a final, desired state, by the gradual alteration 
of the system Hamiltonian. The algorithm requires a Hamiltonian (^^a^ 

whose ground state is easily found, and a Hamiltonian (^Hp^ whose ground 

state encodes the solution to an appropriate mathematical problem. The 
adiabatic algorithm consists of three steps: prepare the ground state of Ha, 
vary the Hamiltonian according to 



with r sufficiently large so that the process is approximately adiabatic, and 
perform an appropriate measurement on the final state, which encodes the 
solution to the problem. 
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This method can be straightforwardly apphed to energy minimization 
problems, such as protein folding. Perdomo, et al. [15] produced an in- 
stance of the adiabatic algorithm which is designed to find low-energy con- 
figurations for lattice-based hydrophobic-polar protein models. While it is 
unknown whether this model can be solved for arbitrary cases by any com- 
puter, quantum or classical, it is reasonable to conclude that a quantum 
algorithm would be of use, as the process of protein folding is inherently 
quantum-mechanical fi6\ IT7]. 

The hydrophobic-polar model of protein structure [18] is a simplified 
model of protein folding, taking the individual peptides to be either hy- 
drophobic (H), or polar (P). These peptides are distributed in a chain on 
a 2D or 3D lattice, and an energy is assigned to the configuration of the 
chain by counting the number of adjacencies between P peptides. This en- 
ergy assignment based on spatial co-ordinates is mapped to a Hamiltonian 
consisting of three terms; an on-site term which penalizes conformations 
that have multiple peptides on the same site, a pairwise interaction term 
that rewards assignments with adjacent hydrophobic peptides, and a struc- 
ture constraint Hamiltonian that penalizes spatial conformations that do not 
form chains. Since the positions of the peptides are the quantities of in- 
terest, a measurement in the computational basis is sufficient; no advanced 
measurement technique is required. 

We have seen that quantum information processing methods can be of 
great benefit to theoretical chemistry, providing properties of ground states, 
as well as reaction parameters. In the next section, we describe in further 
detail the requirements for a quantum computer. This is followed by a report 
on the current state of QIP implementations using nuclear and electronic 
spins in an external magnetic field. These NMR and ESR experimental 
approaches are well-suited to digital QIP, but not adiabatic methods, so we 
will focus on digital methods in the remaining sections. 

1.3 NMR and the DiVincenzo criteria 

In order to perform one of the simulation algorithms mentioned above, there 
must first be a system which displays the unique properties of quantum me- 
chanics, while allowing experimental control and measurement. This idea 
has been formalized into five well-accepted requirements known as the Di- 
Vincenzo criteria [IH]: (1) a scalable physical system with well-characterized 
qubits, (2) the ability to initialize the register to a simple fiducial state, such 
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as |0) (3) a universal set of quantum gates, (4) a qubit-specific measure- 
ment capability, and (5) decoherence times much longer than the average gate 
operation time. These criteria are partially satisfied by NMR QIP, where al- 
gorithms are implemented on nuclear spins in molecules, subjected to a large, 
constant external magnetic field and radio-frequency (RF) pulses [20]. This 
system has the advantage of being well-tested arena for the development of 
quantum control techniques in the few-qubit regime (up to 12 qubits as of 
this writing |21j). Below, we describe the relationship between the desired 
properties of a quantum information processor and the particular features of 
the NMR system. 

1.3.1 Scalability with Well-Characterized Qubits: 
Spin-1/2 Nuclei 

In order for a quantum information scheme (an algorithm or a processor) to 
be scalable, the resources it requires must grow only polynomially with the 
size of the system. Any implementation of QIP using k two-level subsystems 
(referred to as qubits) is scalable in the total energy of the system, or equiv- 
alently, the precision with which the energy is measured. NMR QIP employs 
such an architecture, as each spin-1/2 nucleus has two well-defined energy 
levels in a magnetic field, described by a two-dimensional Hamiltonian, 



where Bq is the external magnetic field, 7 is the nuclear gyromagnetic ratio, 
and 6 is the chemical shift term, imposed on each nucleus by its local molec- 
ular environment. When nuclei have distinct 6, they provide qubits that can 
easily be individually addressed through RF pulses. If this is not the case, 
more advanced techniques are required to address the qubits. 

NMR spectroscopy has been widely used in analytic chemistry and other 
disciplines for decades before quantum information research began, result- 
ing in a large body of literature dedicated to measuring and modelling the 
NMR spectra for a given molecule [22l |23]. The techniques developed in 
conventional NMR can be used to determine the spin Hamiltonian, and thus 
characterize qubits, to a sufficient level of precision for quantum informa- 
tion. However, the number of qubits in NMR QIP systems is constrained in 
a practical sense, because the chemical shifts do not grow with the number 
of qubits, limiting the available frequency space for qubit addressing. Using 



H = -{hjBo + 6)Z = -{hjBo + 6) 



1 

-: 
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solid-state NMR techniques, it is possible to use a magnetic field gradient to 
render the qubits distinguishable without the need for distinct 6. However, it 
is difficult in practice to achieve gradients large enough to provide frequency 
resolution for nearby spins. 

Liquid-state NMR experiments have successfully demonstrated universal 
control over NMR systems ranging from a few qubits in the late 1990s and 
early 2000s [25] to up to twelve qubits more recently [21]. While this 
number is not expected to appreciate greatly, NMR QIP experiments with 
up to ~ 20 qubits are likely achievable. The few-qubit experiments we have 
today are sufficient for testing small quantum simulations and general con- 
cepts of QIP. Perhaps more importantly, the techniques developed in liquid- 
and solid-state NMR QIP are exportable to other, more scalable quantum 
information architectures, on which simulations of the type discussed in this 
chapter that go beyond the capability of the best classical computers could 
be performed. 



1.3.2 Initialization: The Pseudopure State 

To perform a quantum algorithm, we need the register to begin in a known 
state such as |0)®". In NMR, the Boltzmann distribution for a single spin has 
a bias toward the ground state of a ~ tanh (^T-^o/fc^r). In ideal conditions, 
it would be possible to extract a pure state in the liquid state from a multi- 
qubit Boltzmann distribution. However, this is difficult in practice. Instead, 
it is possible to produce a multi-qubit pseudopure initial state, 

Ppscudopurc =^/ + «(|0)(0|)^r (7) 

For typical values of 7 and Bq, a is on the order 10~^. Unitary operations 
due to applied pulses will alter the (|0)(0|)'^"' term just as they would a pure 
state, while leaving the identity term unchanged (neither is the identity term 
measurable, as it does not contribute to the total spin magnetization). As 
long as a is large enough to produce a measurable NMR signal, the final 
pseudopure state can be measured. 

As the size of the system increases, the methods that have been used 
to produce pseudopure initial states will result in values of a that decay 
exponentially in the system size, precluding scalability in liquid state NMR 
QIP. However, experimental procedures such as dynamic nuclear polarization 



and algorithmic cooling, to be discussed in Section 2.3, could increase nuclear 
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polarizations to values near unity in suitable solid state systems that include 
electron spins. 

1.3.3 A Universal Set of Quantum Gates: RF Pulses and Spin 
Coupling 

It is possible to implement an arbitrary n-qubit unitary gate using only one- 
and two-qubit gates [26j. Single-qubit gates are realized in NMR QIP by ex- 
erting a time-dependent, radiofrequency (RF) magnetic field over the sample 
in the X — Y plane (the Z direction being determined by the background 
static magnetic field), which implements the rotation operators Ux{a) and 
Uy{P) (or a rotation along any axis in the X — Y plane of the Bloch sphere) , 
where a and /3 are arbitrary angles determined by the strength and duration 
of the RF pulses. Placing unitaries of this type in sequence can produce 
any possible single-qubit operation. Conditional logic (also known as if-then 
logic) in NMR QIP is achieved by allowing the nuclear magnetic state to 
evolve under coupling mechanisms present in the molecule. These coupling 
mechanisms are direct dipolar coupling and electron-mediated dipolar cou- 
pling (i.e. indirect dipolar coupling ), the latter also known as J-couphng. 

One can also express the universality of an implementation of QIP directly 
in terms of the Hamiltonian, in contrast to examining the unitary operators 
generated by the Hamiltonian. This is more useful for quantum simulations, 
since the class of Hamiltonians which can be simulated by a physical system 
is easily obtained from the form of the natural and control Hamiltonians ■ 
Universal control over a system, in this sense, is the ability to simulate any 
Hamiltonian with the same number of degrees of freedom. Importantly, any 
coupled Hamiltonian, together with universal single-qubit control, can be 
used to simulate any other coupled Hamiltonian that has the same degree of 
connectivity ^7] . 

Since all of the spin-^ nuclei in a molecule possess magnetic moments, 
they will possess a pairwise direct dipolar coupling term in the NMR Hamil- 
tonian: 

£tDD /^o^jTfc/o/^ - - \ - -\ fQ\ 

Hjk (3(^i ■ ejk){(Tk ■ ejk) - (Tj ■ at), (8) 

jk 

where rjk is the distance between the two nuclei, a = [X, Y, Z] and ejk is a 
unit vector along the line connecting the two nuclei. In the liquid state, Qjk 
(the angle between ej^ and the external 5-field axis) is, to a good approx- 
imation, uniformly distributed by fast molecular tumbling. This results in 
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zero average interaction strength. Therefore, direct dipolar couphng is only 
useful in systems in which the Qjk are not uniformly distributed, such as solid 
state NMR (where all of the 6^^ are fixed by crystal orientation) and liquid 
crystal NMR (where the Qjk are distributed non-uniformly) . In addition to 
the direct dipolar coupling, there exists an indirect coupling, known as the 
J-coupling. This is a relatively weak coupling mediated by the electron cloud 
in the molecule, according to 

Hi = 2na, ■ J,k ■ al (9) 

where Jjk is a 3 x 3 real matrix. This coupling term does not vanish, even 
under rapid tumbling. In isotropic liquids, it has the Heisenberg form aj ■ ak, 
which reduces to ZjZk when the two spins have a chemical shift difference 
much larger than the value of J; this is the typical generator of a two-qubit 
gate in liquid-state NMR. The variety of natural coupling terms renders QIP 
possible on many different molecules, in the liquid, solid, or liquid crystalline 
states. 



1.3.4 Measurement: Free Induction Decay 

The apparatus (an RF coil) used to implement rotations about axes in the 
X — Y plane can also be used to detect ensemble magnetization signals from 
the sample in the X — Y plane. In order to translate the logical state of an 
NMR QIP register into an X observable value, a readout pulse of angle | is 
applied about the X-axis at each qubit resonance frequency. This results in 



a signal, proportional to the pseudopure parameter a (see Subsection 1.3.2), 
from which the logical state can be extracted. This process does not result 
in projective measurement, but still allows universal computation including 
full state tomography [28]. In addition, the expectation value of any Hermi- 
tian operator can be obtained in a single shot, allowing the measurement of 
arbitrary operators by quadrature detection 



1.3.5 Noise and decoherence: 

T1/T2 Vs. Coupling Strength 

In order to perform an algorithm, the time used to implement the appropri- 
ate pulses must be much shorter than the characteristic timescale of deco- 
herence. In addition to the noise related to imperfect implementation of the 
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Ci(Hz) 


C2(Hz) 




T*(ms) 


Ti(s) 




5693 


237 


828 


2.4 


162 


C2 




1748 


1020 


2.6 


326 


c 






-3358 


3.1 


314 



Table 1: (from [22]) Characteristic timescales for control in solid-state 
NMR QIP, using the ^^C-labelled malonic acid molecule, can be gleaned 
from the dipolar coupling strengths (off-diagonal elements in second and 
third columns) and chemical shifts (diagonal elements), and compared with 
relaxation and dephasing times (rightmost columns). Note the description 
of decoherence in terms of , which combines decoherence from inhomoge- 
neous and homogeneous sources. Much longer T2 ~ 100 ms was observed 
with suitable dipolar refocusing / dynamical decoupling sequences applied 



gates (these are typically 'coherent' errors), the unavoidable interaction with 
the surrounding environment leads to true decoherence. In NMR, there are 
two main mechanisms for this decoherence: thermal equilibration (energy re- 
laxation) and dephasing, characterized by timescales Ti and T2*, respectively. 
T2* represents the ensemble, inhomogeneous dephasing time of the system, 
whereas T2 is the intrinsic decoherence time of each single quantum subsys- 
tem. T2* < T2 can often be improved to T2 by applying decoupling pulses as 



described in Section |2.2[ In NMR QIP, the coupling terms between qubits 
are usually the smallest terms in the Hamiltonian, leading to relatively slow 
two-qubit gates. The solid-state malonic acid system provides an example, 
shown in Table [Tj the carbon-carbon coupling ranges from 200 Hz to 1 kHz. 
This indicates that the associated controUed-Z gate requires a time of 1/2 J ~ 
0.5 - 2.5 ms, allowing the implementation of several two-qubit gates before a 
single decay period has elapsed. This suffices for simple algorithms, but 
a larger ratio of decoherence-to-gate time is required in order to implement 
more complex algorithms, and in particular, to implement error correction. 
If the noise is below a threshold [HT], it is possible to use fault-tolerant meth- 
ods to counteract its effects. If the noise exceeds this threshold, control must 
be improved until the threshold is attained. It is possible to assess the noise 
level using benchmarking [321 [33]. In an NMR system where control was 
optimized, it has been possible to reach an error per gate rate of 10~^ for 
single-qubit gates and 10~^ for two-qubit gates 
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Ongoing improvements to the implementation of one- and two-qubit gates 
in magnetic resonance QIP have arisen, in part, due to: 

• optimization of RF pulses to produce high-fidelity unitary operations, 

• isolating the system of interest from unwanted degrees of freedom, and 

• taking advantage of the unique properties of the electron spin. 

In the following sections, we describe recent progress in these directions. 

2 Quantum Control in Magnetic Resonance 
QIP 

2.1 Advances in Pulse Engineering 

In NMR QIP, pulse engineering is the practice of developing control tech- 
niques for either generating coherent transfer from an initial spin state to a 
desired spin state (state-to-statc transfer) or producing an effective Hamil- 
tonian that implements a desired unitary gate U (unitary propagator) by 
manipulating the external RF field. This can be achieved by augmenting the 
internal Hamiltonian with a control Hamiltonian (in the laboratory frame), 

Hc{t) = ^^[cos(c^RFi + (t>k{t))X + sin(c^RFi + (t>k{t))Yl 

k 

where a;fe(t) and 0fe(t) denote the amplitude and the phase applied at a 
transmitter frequency WRp-The time evolution is typically divided into a set 
of steps. At each step j, the evolution is given by the unitary propagator 
Uj = exp{—iAtHj{t)}, where Hj{t) is the total Hamiltonian with a control 
set {oul, (pi}. The final density operator is given by p(T) = (Jm---UiPoUi---U'Ii- 
The task of pulse engineering is to find a set of {u;^, 0^} so that the resulting 
dynamics are sufficiently close to those desired. The quality of the pulse is 
evaluated by a fidelity function $ that is proportional to overlap between the 
obtained operator and the target operator. 

Pulse engineering techniques must address experimental hmitations for 
practical applications: qubit selectivity due to finite frequency bandwidth, 
instrumental errors (such as static and RF field inhomogeneity), miscalibra- 
tion of pulse power or duration, and frequency offset. 
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Pulse design must take these artefacts into account in order to achieve 
robust pulses with high fidelity on the ensemble. This causes a substantial 
increase in the complexity of the pulse design problem. For example, errors 
caused by chemical shift dispersion, RF inhomogeneity and RF miscalibration 
can be suppressed if the parameters are sampled over a range of discrete 
values determined by the uncertainty. Then the total fidelity function to be 

optimized can be measured as $tot = ^^^(a^a), where X (y^ IS 9j particular 

a 

value of some parameter seen by a fraction of spins, e.g. RF field amplitude. 
In this case, the optimization problem becomes significantly more difficult, 
since the effective parameter space is much larger. 

Over the last 30 years, numerous techniques have been developed for 
control in NMR QIP. Traditionally, average Hamiltonian theory (AHT) has 
been a powerful tool that provides intuitive guidelines for constructing pulse 
sequences in relatively simple cases [35] . In particular, composite pulses [36] , 
adiabatic pulses [371 |3H1 EH 00] , and shaped pulses [H] were introduced dur- 
ing the earlier development of NMR QIP to better compensate for static and 
RF field inhomogeneity by increasing the number of degrees of freedom in 
the pulse shape. However, the long pulse times produced by these techniques 
lead to greater decoherence and relaxation effects, and interference of selec- 
tive pulses simultaneously applied to different spins. However, not all errors 
can be corrected using these techniques [12] . Strongly modulated pulses [12] 
average out unwanted evolution by using strong control fields to drive the 
spins, so that a desired unitary operator can be synthesized directly. This 
method uses high-power pulses that decrease the required pulse duration, 
and hence reduce the effect of decoherence and relaxation. Penalty func- 
tions on power, frequency and time period can be applied to obtain more 
experimentally acceptable solutions. Moreover, it allows for the incorpora- 
tion of robustness against errors caused by static or RF field inhomogeneities. 
Strongly modulated pulses with ideal fidelities of order 99% are performed 
for universal control of liquid-state NMR system with up to six qubits [12] . 
The method was also studied in a three-qubit solid-state NMR system, single 
crystal malonic acid. Universal control with ideal fidelities 98% was achieved 
in the strong-coupling regime of the direct nuclear dipole-dipole couplings, 
in which the state evolutions under the internal Hamiltonian are more com- 
plex than those of typical liquid-state NMR system [30j. In addition, the 
transformations generated by the strongly modulated pulses were shown to 
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be robust against the ensemble inhomogeneities that are present in sohd- 
state NMR systems [30]. The cost of this method (as well as the optimal 
control methods described below) is that the full system dynamics must be 
simulated in the course of finding pulses, so that even using clever methods 
such as gradient-based approaches to improve efficiency does not render the 
method scalable to arbitrary numbers of qubits. 

Recently, optimal control theory (OCT), originally developed for prob- 
lems in engineering [ISj HH |l5|, has been used for systematic optimization 
of pulse designs in NMR QIP. Analytical solutions to time-optimal real- 
ization of unitary operations can be obtained by formulating a variational 
principle to reduce the problem to a set of first-order ordinary differential 
equations [IHl H?]. For more than two qubits, the analytical approach is 
intractable, and numerical methods are required. 

For larger systems, the majority of control methods for optimizing the 
efficiency of coherent transfer are based on a gradient approach, such as the 
Gradient Ascent Pulse Engineering (GRAPE) algorithm [IB]. The GRAPE 
algorithm can be summarized as follows. It begins with an initial guess for 
a set of control parameters Uk{j) for all time steps j < N. Then in each 
iteration the algorithm evaluates S^tot/Suk{j) for j < N and updates Uk{j) 
as: 

Uk{j) ^Uk{j) + e- — --. 

The iteration continues until the improvement in the performance index ^tot 
is smaller than a chosen threshold value. GRAPE also evaluates penalty 
functions to reflect realistic experimental constraints. The gradient approach 
enables GRAPE algorithm to calculate full time evolutions much faster than 
conventional (brute force) numerical difference methods. Consequently, the 
number of pulse parameters to be optimized can be orders of magnitude 
larger than in conventional approaches. A number of recent NMR QIP ex- 
periments [m EOl EH [5211531 IMl ES] have utilized GRAPE pulses to achieve 
high control fldelities. However, just as with any optimal control method, the 
GRAPE algorithm cannot deterministically locate global minima but instead 
flnds local minima in the search space. Therefore the ultimate fldelities that 
can be achieved are limited by the initial guess, even if unit fldelity pulses 
are possible in principle. Figure [T] illustrates an example of a path taken by 
the spin magnetization vector on the Bloch sphere during a GRAPE pulse 
designed for a | rotation around the X-axis. 

Another well-developed numerical method for quantum control is Krotov- 
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Figure 1: Trajectory of the spin magnetization vector on the Bloch sphere 
for a pulse found by the GRAPE method [IS] . The pulse is designed for a | 
rotation around the X-axis. 

based numerical method [561 E3 EHl EH] . The Krotov-based method allows 
large changes in control parameters from one iteration to the next, and im- 
mediately exploits all available information at each time step, monotonically 
improving the objective fidelity functional at each iteration. Maximov et 
al. [IS] analysed the Krotov-based algorithm in the context of NMR spec- 
troscopy and compared it to GRAPE, concluding that the Krotov-based 
algorithm consumes less computational resources per iteration and is much 
better than GRAPE for making an initial guess towards a global maximum. 
However, the larger step size in the Krotov-based method limits efficiency as 
the solution gets closer to the desired fidelity. There is an open question as 
to whether high efficiency can be obtained by combining the two methods, 
using the Krotov algorithm to quickly prepare a good initial pulse sequence 
to load into GRAPE for final refinement f?5]. 

As mentioned above, a major drawback with these methods is that com- 
putations cost grows quickly with increased system size. As quantum pro- 
cessors become larger, new techniques will need to be developed that might 
allow these OCT methods to be applied within blocks of nearby qubits (in 
frequency space) and scalable pulse design techniques between blocks [60j . 

2.2 Advances in Dynamical Decoupling 

In order to achieve high-fidelity control in NMR QIP, one must address 
system-environment interaction, as well as unintended evolution under in- 



16 



ternal couplings between spins. The coupling interaction is "always on" and 
therefore must be suppressed to perform single qubit gates. This process 
is called refocusing and is relatively simple in liquid state NMR, since the 
only coupling term is of the form Z^Z]^. Solid-state NMR Hamiltonians have 
more complex coupling terms, and are therefore more difficult to refocus. 
Techniques based on average Hamiltonian theory, such as magic angle spin- 
ning |6T], Lee-Goldburg decoupling [62] and multiple-pulse techniques [35] , 
have proven useful in refocusing these terms. Furthermore, recent devel- 
opments in optimal control theory (see Section 2.1) allow the refocusing of 
unwanted internal interactions for arbitrary coupling forms, as long as the 
Hamiltonian is well-known. 

Hereinafter, we narrow the discussion to system-environment, or dynam- 
ical decoupling (DD). A DD scheme is a sequence of control fields applied to 
a system for some time with the objective of increasing coherence time by 
attenuating the system-environment interaction. Starting from Hahn's dis- 
covery of spin-echo in 1950 [03] , many methods to accomplish this have been 
actively studied, one seminal example being the Carr-Purcell-Merboom-Gill 
(CPMG) sequence [M] which is successful at suppressing single-axis, low- 
frequency noise. The performance of realistic DD sequences is limited by 
experimental imperfections (see Section 2.1), so it is important to design a 
DD scheme that is more robust to the most relevant errors. 

In traditional DD, a pulse sequence is periodically applied to reduce un- 
desired terms of the system-bath interaction Hamiltonian, known as periodic 
dynamical decoupling (PDD). These sequences were improved in 2005 by 
Khodjasteh and Lidar [65], who introduced concatenated DD pulse sequences 
(CDD), such as = PnXpnZpnXp^Z , where Pn = tq is a delay between 
pulses and n is a concatenation level. CDD is significantly more efficient 
at decoupling than PDD with equal pulse numbers and is more robust to 
both random and systematic control errors. CDD can reduce the system- 
environment interaction to order where t is total duration of the cycle 
and n is the concatenation level. The caveat of this method is that due to the 
nature of concatenation the number of necessary pulses grows exponentially 
(4") in the concatenation level, whereas order of suppressed error only grows 
linearly. 

In 2007, Uhrig introduced an optimized dynamical decoupling sequence 
based on gradient moment nulling in NMR [66], called Uhrig Dynamical 
Decoupling (UDD) ESI EH]- UDD was first introduced for a specific 
spin-boson model, and was later discovered |68l [69] to be system indepen- 
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dent. What distinguishes UDD from conventional DD schemes is that the 
time delay blocks between n pulses are not equal. Uhrig showed that split- 
ting the total time interval t into smaller intervals 0, 6it, 5„t (where Sj = 
sin^[7rj/(2n + 2)]), can suppress decoherence up to order without expo- 
nential cost. Moreover, Uhrig's simulation showed that the performance of 
decoupling from the environment becomes independent of system-bath cou- 
pling strengths for long sequences. One drawback of UDD is that it works 
for only single-axis error, so for example, cannot simultaneously suppress Ti 
and T2 processes. 

In an attempt to refocus the effects of noise about multiple axes, Uhrig 
concatenated the UDD sequence (CUDD) [TO]. The total number of pulses 
to suppress decoherence and relaxation to order is proportional to 2", 
decreasing the resources by 2^ from CDD sequences. West et al. further 
improved on this by creating a new UDD-based sequence that suppresses 
decoherence and relaxation to order n using {n + 1)^ pulse intervals [71] . 
Their scheme integrates two sequences, Ti-correcting UDD and T2-correcting 
UDD, and is known as quadratic UDD (QDD). 

UDD-based methods are optimal when the noise has a sharp high-frequency 
cutoff [72l [731 ES [ZH [75] . However, in a low-frequency-dominated noise en- 
vironment, conventional CP (or CPMG) is preferred [73l[71l[75]. Borneman 
et al. [TB] pointed out that, although the CPMG sequence is inherently tol- 
erant of field inhomogeneities and pulse calibration errors, the robustness of 
the sequence is limited by the quality of the RF pulses used. They adapted 
the GRAPE algorithm (see Section 2.1) and applied iterative optimization 
method to design general refocusing pulses (universal vr rotation around the 
F-axis) with high fidelity over a wide range of resonance offset frequencies 
and RF amplitudes. Borneman et al were able to find a pulse which refo- 
cuses 99% of the initial magnetization over a range of frequency offsets of 
±10 KHz (four times greater than the maximum RF amplitude) for uniform 
RF field. For RF inhomogeneity of ±10%, they were able to find a pulse 
that refocuses 98% of the initial magnetization over a frequency range of 
±8 KHz (3.2 times the maximum RF amplitude), an improvement over any 
previously published refocusing pulses of similar duration and maximum RF 
amplitude when applied in a CPMG sequence. 

Souza et al. [77] compared the performance of a number of pulse sequences 
in order to find a decoupling scheme that is robust against pulse imperfections 
or control errors. They considered two possible approaches to compensate the 
imperfections: use composite pulses that the error correction is done "inside" 
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the pulse (robust pulse), or create the sequence where the error from one 
pulse is averaged out in subsequent pulses (self-correcting sequence). In the 
experiment a standard pulse (non-robust) and the Knill pulse [78] 7r^/6+</> — 
TT,^ — 7rjr/2+<^ — TT^ — 7r^/6+0, which is robust against flip-angle and off-resonance 
errors are inserted into the CPMG, periodic XY-4 [79], and CDD. Then 
decay times of the magnetization as a function of the duty cycle, the total 
pulsing time in the sequence divided by total duration of the sequence, are 
analysed and compared for each case. A new sequence called Knill Dynamical 
Decoupling (KDD) which is achieved by inserting delays between the vr pulses 
of the Knill pulse, is also examined. For small duty cycles, sequences with 
non-robust pulses are superior, due to the shorter cycle time of the non- 
robust pulse sequence when constant duty cycles are compared. However, the 
performance of these sequences saturates or decreases with increasing duty 
cycle. Comparatively, robust pulses continue to increase in performance for 
large duty cycles. KDD is found to have the best performance for large duty 
cycles, and comparable to sequences without robust pulses for small duty 
cycles. 

2.3 Control in the Electron-Nuclear System 

Many of the techniques used to control spin-y2 nuclei in NMR can also be 
applied to electrons, comprising electron spin resonance (ESR). As the 7 of 
an electron is ~ 660 times greater than that of a proton, this leads to a 
much larger Larmor frequency and higher polarization, but also faster deco- 
herence and relaxation than nuclear states. The coupling between electron 
and nuclear spins is governed by the hyperfine interaction: 

HM = aE-A- a]^, (10) 

where E denotes the electron, denotes the nucleus, and A denotes the 
hyperfine coupling tensor. The combined electron-nuclear solid state spin 
system is appealing for QIP because it is possible to exploit the strengths of 
each type of spin: electron spins possess higher polarization for initialization 
and readout, and can be manipulated on faster timescales, while nuclear 
spins are ideal for long-term storage of coherent states [80J. 

Although ESR provides much better signal-to-noise ratio than NMR, it 
usually suffers from broader linewidths. For an electron-nuclear spin system, 
one is required to characterize its internal Hamiltonian as precisely as pos- 
sible in order to achieve high fidelity control. However, the nuclear Zeeman 
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term in the Hamiltonian is small relative to terms involving the electron, 
and often it is not possible to measure directly in an ESR experiment due to 
spin selection rules ^Sl]- Moreover, the ESR linewidth can be comparable to 
the hyperfine splitting which in turn makes it difficult to study the hyper- 
fine interaction term. Fortunately, the well-established magnetic resonance 
techniques of electron spin echo envelope modulation (ESEEM) [SH [H2] and 
electron nuclear double resonance (ENDOR) spectroscopy [HH ESJ [H] can 
precisely quantify nuclear transition frequencies and hyperfine couplings. In 
ENDOR, both NMR and ESR transitions are driven directly. Thus EN- 
DOR spectroscopy is a promising technique for controlling electron-nuclear 
systems for QIP. However, it should be noted that ENDOR experiments re- 
quire additional RF hardware, and nuclear spin fiip times are limited by the 
small nuclear gyromagnetic ratio and the strength of the RF field that can 
be achieved. 

The anisotropy of hyperfine interaction allows controlling nuclear spins 
solely by irradiating electron spin transitions [SS] . The indirect control tech- 
nique is advantageous: it simplifies the instrument design as an additional RF 
interface is not needed, and faster control of a nuclear spin is achievable when 
the hyperfine coupling strength exceeds the Larmor frequency of the nucleus 
at given external field [86]. The following section reviews this approach and 
recent experiments demonstrating the universal control in electron-nuclear 
systems. 



2.3.1 Indirect Control via the Anisotropic Hyperfine Interaction 

If the coupling frequency due to hyperfine interaction is larger than the nu- 
clear Larmor frequency, efficient nuclear control can be obtained by manip- 
ulating the electron [86], producing an effective magnetic fields seen by the 
nucleus that depends on the electron spin state. With the anisotropic cou- 
pling, the nucleus sees an effective field 

5eff = (Bo ± ^) Z ± -^X, (11) 

where A is the 2;z-component of the hyperfine tensor [81] and B is ^ Al^ + Aly 
the component of the hyperfine tensor in the x-direction on the nucleus, in a 
specially chosen frame. Bq is the static magnetic field and 77V is the nuclear 
gyromagnetic ratio. The ± sign takes the value + when the electron spin is 
parallel to the external field, and — when it is anti-parallel. Any rotation 
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of the Bloch sphere can be generated by repeated rotations about these two 
distinct axes. Universal single-qubit control of the nuclear spin can then 
be achieved indirectly, since the free precession of the nucleus about these 
axes does not require RF pulses resonant with nuclear transitions. Hodges 



et al. implemented GRAPE algorithm (see Section 2.1) and showed univer- 
sal control of a le-ln system (malonic acid radical) [B5]- In [SD], Zhang et 
al. also utilized the GRAPE method to perform an entangling gate between 
two nuclear spins in a le-2n system (singly ^^C-labelled malonic acid radi- 
cal). Schematics of the malonic acid radical, an energy level diagram showing 
creation and detection of double-nuclear coherence, and the double-nuclear 
coherence echo signal from [5U] are shown in Figure [2] These experiments 
represented the first steps toward reaching high-fidelity universal control of a 
one-electron, N-nuclear spin hyperfine coupled system by using the electron 
as an actuator. 

In the following sections, we discuss additional quantum information tech- 
niques which exploit the properties of combined NMR/ESR systems: 

• dynamic nuclear polarization and algorithmic cooling, which increase 
the polarizations of nuclear spins, 

• spin-bus implementations and parallel information transfer, which pro- 
duce two-qubit gates between nuclear spins mediated by electrons. 

These two techniques, described below, reinforce the potential of combining 
nuclear and electron spin control to better satisfy the DiVincenzo criteria in 
NMR QIP. 



2.3.2 Dynamic Nuclear Polarization and Algorithmic Cooling 

A nuclear spin in a magnetic field at room temperature has an equilibrium 
density operator described by a Boltzmann distribution, with a typical bias 
of ~ 10~^. Although this is sufficient to demonstrate the principles of QIP, 
it does not allow for resetting ancilla qubits as needed, for example, in quan- 
tum error correction. Since the nuclear magnetic moment is so weak, this 
problem cannot be simply resolved by brute-force cooling or increasing the 
magnetic field. For example, reducing the temperature to 1 K and using the 
largest available magnetic fields increases the polarization only to the order 
10~^. Therefore, it is necessary to develop novel techniques to provide initial 
polarization that approaches unity. 
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Figure 2: From pQ]. (a) upper: Schematic of the singly ^^C-labelled malonic 
acid radical lower: Table of Hamiltonian parameters for the crystal orienta- 
tion used in the experiment. Larmor frequencies are along the diagonal, and 
hyperfine coupling coefficients are off-diagonal, with all frequencies in MHz. 
Ti and T2 are shown in the two right-most columns in microseconds, (b) 
Schematic of the experimental sequence for creating and detecting a double- 
nuclear coherence in the electron spin down manifold, (c) The double-nuclear 
coherence echo signal versus delay time in experiment (red crosses, solid line) 
and simulation (blue dashed line) for the transition between \l Ih^c) and 
\l Oh^c) in (b) (electron spin-down manifold). 

The much larger electron magnetic moment can be used as a polariza- 
tion source. A swap gate between a nuclear spin and an electron, each at 
thermal equilibrium, acts as a polarization swap and gives the nuclear spin 
a large non-equilibrium polarization. This type of process is called dynamic 
nuclear polarization, or DNP. In the context of QIP, such polarization trans- 
fers have been implemented in ^^NtQCgo ^7], achieving a single-spin nuclear 
polarization of 62% at 4.2 K, in an 8.6 T field. Also, DNP has been used in 
conjunction with spin diffusion in silicon microparticles to produce bulk spin 
polarizations up to 5% at 1.5 K and 2.35 T p8]. Nuclear Ti for these systems 
was shown to be dependent on the particle size, making the polarization de- 
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cay controllable in principle. The ^^N nucleus at a Nitrogen vacancy (NV) 
centre has been initialized to a polarization of 98% by optically pumping the 
NV defect [89], and similar polarizations can be reached for proximate ^^C 
nuclei. These examples show that control of a single electron spin coupled 
to an NMR QIP register can be used to great advantage. 

It is possible to obtain nuclear qubits with still greater polarization by an 
implementation-independent method, algorithmic cooling. This technique 
combines reversible entropy compression (unitary process) and open quan- 
tum systems cooling (non-unitary process). For an m-qubit system with each 
qubit possessing equal entropy, a reversible entropy compression process can 
compress entropy into m — n qubits, increasing their effective spin temper- 
ature while cooling the subset of n qubits. Here, lower entropy corresponds 
to higher polarization. The system is in contact with a large heat reservoir 
(bath) whose entropy is lower than the initial entropy of each system qubit. 
SWAP gates are then applied to exchange the system entropy and the heat 
bath entropy to cool down m — n qubits that were heated during the com- 
pression step. The local part of the heat bath quickly relaxes back to its 
equilibrium state and reacquires its initial high polarization. This two-step 
process is applied iteratively until no further cooling is possible. Given m 
qubits and a heat bath polarization of e^, > 2"™, it is possible to cool the 
system very close to its ground state with only polynomially increasing re- 
sources [9D]. However, for < 2""^, the maximum polarization attainable 
by a single spin is eb2"*~^ [29]. Algorithmic cooling has recently been imple- 
mented in liquid and solid state NMR [911 [501 192] . Since electron spins possess 
higher polarization and much have shorter Ti times than nuclear spins at a 
given magnetic field strength and temperature, having nuclear qubits able 
to interact with an electron spin 'bath' is a promising path to nuclear qubit 
initialization. Single quantum systems, such as NV centres with a sufficient 
number of coupled ^^C spins, would be ideal for demonstrating nuclear qubit 
purification by algorithmic cooling. 

2.3.3 Spin Buses and Parallel Information Transfer 

Since a single electron can couple to multiple nuclei, it can be used to trans- 
fer information between them, creating an effective coupling. This indirect 
coupling is the basis of the S-bus [93], useful for performing multi-qubit 
gates when the electron-nuclear and electron-electron coupling are much 
stronger than the nuclear-nuclear coupling. This concept was used to per- 
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form Deutsch's Algorithm [^2], using a system of two nuclear qubits coupled 
by an electron in CaF2:Ce^+. 

The indirect coupling of nuclei using information transfer of the type 
described above has two flaws which must be overcome to confirm its utility 
as a nuclear control method. The first is that the state being transferred is 
subject to strong decoherence and relaxation while it is stored on the electron. 
The second is that only two nuclei can be coupled through the bus at the 
same time; two-qubit gates cannot be performed in parallel. These problems 
were recently resolved [l94J, where the assumed architecture consists of two 
local nodes (taken to be sets of n nuclei, each coupled to an electron via 
anisotropic hyperfine interaction), with the only inter-node coupling (dipolar 
or exchange interaction) being between electrons. By using the interaction 
frame, Borneman, Granade and Cory showed that states of the multi-nucleus 
nodes can be swapped in parallel, effectively performing n two-qubit gates 
simultaneously. Also, the effect of decoherence on the electron is mitigated 
by ensuring that no computational state is stored on the electron. Figure [3] 
illustrates the idea presented in |Mj . 




Nodel Node 2 



Figure 3: Schematic of 2 x (le-3n) Node: The nodes are taken to be identi- 
cal, with resolved anisotropic hyperfine interactions (solid red lines) between 
electron actuator spins and nuclear processor spins. The local processors are 
initially disjoint, but may be effectively coupled (dotted lines) by modulating 
an isotropic actuator exchange or dipolar interaction (solid blue double line) 
and moving into an appropriate microwave Hamiltonian interaction frame. 
The spin labelling is for electron actuator spins and Uij for nuclear proces- 
sor spins, where i and j label the nodes and the qubits, respectively. |9l] 

The electron spin, when used as a component of the NMR QIP sys- 
tem, possesses properties that complement those of the nuclear spins. While 
the nuclear spin has a longer coherence time and low initial polarization. 
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the electron has a short coherence time and high initial polarization. The 
electron-nuclear coupling also permits an array of techniques, making the 
electron a valuable asset to NMR QIP. 

An example of electron-nuclear system application for QIP that has been 
extensively investigated recently is entanglement generation [OSl ESI EZl EHl 
inSl 1100] - Preparing entangled states in NMR QIP is impracticable since 
required spin polarization according to the positive partial transpose (PPT) 
criterion |10H I102j for entanglement is far above what is reachable. Electron- 
nuclear QIP at high field and low temperature, in conjunction with optimal 
control techniques, can potentially provide a highly entangled quantum state 
that is a key ingredient in many quantum algorithms and applications. 

3 NMR QIP For Chemistry 

3.1 Recent Experiments in NMR Quantum Simulation 

The advances described in the previous section can enable the development 
of larger, more precise processors for quantum information. We focus now 
on three recent experiments that display the current capabilities of the NMR 
QIP system, and prove its utility as a testbed for future implementations. 
These experiments are motivated by the study of physical systems rather 
than chemical processes, however, the methods within are easily adapted to 
problems of chemical interest. 

3.1.1 Simulation of Burgers' Equation 

One challenge which presents itself to scientists of many disciplines is the 
solution of non-linear differential equations. Numerical solution of these 
equations is costly, this cost having motivated the development of numer- 
ous simulation methods. One such method, which takes advantage of quan- 
tum resources, is the quantum lattice gas algorithm on a type-II quantum 
computer. 

Type-II quantum computers |104] are networks of small quantum proces- 
sors which communicate classically to form a large analog computer. Oper- 
ations that can be carried out in parallel are well-suited to implementation 
on such a processor, since the individual QIP nodes are distinct and can 
therefore be addressed simultaneously. The quantum lattice gas algorithm 
consists of such operations; it is divided into three steps which are iterated: 
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Figure 4: From |1U3] . Comparison of analytic and numerical solutions to 
Burgers' equation. Numerical data has been obtained on a 16-node Type-II 
quantum computer. 

1. A state is prepared on the individual QIP nodes, corresponding to 
initial (or transient) conditions of the dynamical model being studied. 
The state on a given node is 

(v/l-/a(x,,to)|0) + ^fa{xi,tomr' (12) 

denoting two copies of the state encoding the value of the function to 
be propagated. 

2. A collision operator U is applied in parallel across all of the nodes, 
which simulates a single finite time step for the dynamics being mod- 
elled. 

3. A measurement is made on the nodes and the resulting transient state 
is fed back to step 1. Classical communication is used, in this step, to 
transfer the results of step 2 between QIP nodes. 

The resulting class of finite difference equations which can be numerically 
solved by this method is that of the form 

fa{xi+ea,tn+l) = fa{xi,tn) + {i'{xi,tn)\U'' flaU - na\i'{xi,tn)) (13) 
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where = ±1, U is the colhsion operator, and fa is the function whose 
evolution is being simulated. 

Discretization of Burgers' equation yields a finite difference equation of 
the form described above, it is a non-linear differential equation which de- 
scribes shock formation in fluid dynamics: 

dtu{z, t) + u{z, t)dzu{z, t) = udzzui^z, t) (14) 

Due to the simplicity of this equation, it can be solved analytically. From |103] , 
the analytical solution to Burgers' equation is compared to the simulation 
on a 16-node type-II NMR quantum processor, shown in figure |4| 





Ci 


C2 


C3 


Ci 


M 


H, 


H2 




-1914.06 


40.5 


1.5 


7 


127 


3.9 


6.3 


c. 




-18115.10 


69.9 


1.3 


-7.1 


155.1 


-0.6 


Cs 






-15157.41 


73.2 


6.6 


-1.8 


163 


C4 








-21148.90 


-0.9 


6.5 


3.6 


M 










230.43 


6.9 


-1.7 


Hi 




n 








-2370.80 


15.5 














-1774.47 



Figure 5: From |105] . Spin Hamiltonian parameters and molecular diagram 
for transcrotonic acid. 

3.1.2 Simulation of the Fano- Anderson Model 

Another challenge with which physicists are often presented is that of simu- 
lating the evolution under a Hamiltonian representing the energy landscape 
of a large number of identical particles. The Fano- Anderson model possesses 
such a Hamiltonian, acting on n spin-less fermions constrained to a ring, 
surrounding an impurity: 

n-l 

^ = E '^A^'^^ + ^^^^ + ^(4o& + (15) 
1=0 

where Cfc, is a fermionic annihilation operator on the conduction mode ki, his 
a fermionic annihilation operator on the impurity in the center, c|.^ and are 
the respective creation operators, and e^j , e and V denote the strengths of the 
terms in the Hamiltonian corresponding to occupation of a site, occupation 
of the impurity and tunnelling between sites and the impurity, respectively. 
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This Hamiltonian is exactly diagonalizable, however, its simulation is an 
important step in validating QIP methods for application to simulations in 
condensed matter theory. 
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Figure 6: From |inSj. The kernel of the quantum simulation is the implemen- 
tation of a unitary operator that evolves the qubit states according to the 
desired Hamiltonian. Given the capabihty to produce controlled operations 
between an ancilla qubit and the simulation register, correlation functions of 
interest can be directly measured. 

The simulation performed by Negrevergne, et al. |1U5] concerns the small- 
est possible Fano-Anderson model, one site interacting with an impurity. 
This algorithm requires three qubits, one for the site, one for the impurity 
and one ancilla qubit, required for indirect measurement. It is implemented 
in a liquid-state NMR system, using transcrotonic acid, whose Hamiltonian 
parameters are shown in Figure [5j 

Again, the simulation is divided into three broad steps. However, iteration 
is unnecessary in this experiment, since the Hamiltonian to be simulated can 
be exponentiated continuously. First, the system is prepared in an initial 
state, based on the pseudopure state. Next, in order to perform an indirect 
measurement, a two-qubit gate couples the ancilla to the simulator, which 
is then evolved according to the Fano-Anderson Hamiltonian, and decoupled 
using another two-qubit gate, shown in Figure [6] 

Two quantities of interest were obtained in this experiment. G{t), the 
correlation between the states b'^{t)\FS) and b'^{0)\FS), where \FS) is the 
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state corresponding to the filled Fermi sea. In the NMR QIP system, this 
correlation function is {10\ exp (iHt) a^exp (^—iHt) aUlO). The results of 
the measurement of this function in the simulation are shown in Figure [7) 
Another quantity, of more general interest, is the spectrum of the Hamilto- 
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Figure 7: From |105j . Analytic and numerical solutions for the correlation 
function G{t), whose evaluation is given diagrammatically in Figure [6j 

nian. This can be measured with a simpler quantum circuit, and results in 
excellent agreement with theory were also obtained in reference [105] . 

3.1.3 Simulation of Frustrated Magnetism 

Any spin system which has interaction Hamiltonians that cannot be simulta- 
neously minimized is geometrically frustrated. Frustration is a fundamental 
problem in the study of magnetism, since its presence means that the Hamil- 
tonian cannot be divided into small subsystems which can then be studied 
individually to obtain information regarding a global property. This exac- 
erbates the central problem of spin system simulation, i.e. the exponential 
number of degrees of freedom required to simulate such a system. Zhang et 
al. [r06j have performed a digital quantum simulation of the fundamental 
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Figure 8: From |106j . (a) Spin configurations of tfie Ising Hamiltonian at 
zero temperature and field h. There is a six-fold degeneracy in tlie ground 
state, leading to a non-zero entropy, (b) Tfiermal pliase diagram for tfie same 
Hamiltonian sliowing total magnetization versus field and temperature. 



building block of such a frustrated system: the three-spin anti-ferromagnetic 
Ising model: 

H = J{ZiZ2 + ^2^3 + ZiZ^) + h{Zi + Z2 + Z3), (16) 

where J > 0. This simple Hamiltonian can be solved analytically, due to 
its small size. The central technique in this work is to simulate an arbitrary 
thermal state using pseudopure state preparation. In this manner, the phase 
diagram in h, J and T can be explored, deriving multiple properties of the 
resulting states. The pseudopure portion of the prepared state is 

1^^) = J2 v/exp(-/3Efc)/Z|0,) (17) 

k 

where the are the eigenstates of the Hamiltonian with energy E^, Z is the 
partition function and /3 = 1/T. The two physical quantities which Zhang, 



30 



et al. focus on are the total magnetization of the system Zi + Z2 + Z^, and 
the entropy of the resuhing state S = Tr(p^ log p/3), where p/3 is the thermal 
density matrix of the Ising system (see Figure [9]). A peak in entropy is 
observed at /i = indicating frustration of the magnet, in agreement with 
theory. 

These three experiments serve to demonstrate the validity of the princi- 
ples of quantum simulation, and demonstrate the necessary control to achieve 
scientific results, in principle. The greatest challenge to implementation of 
these algorithms on a larger scale is the inability to control very large Hilbert 
spaces. The advances presented in section 2.3, for example, may be helpful 
in obtaining this control. 



4 Prospects for Engineered Spin-Based QIPs 

Nuclear and electron magnetic resonance on bulk materials containing nat- 
ural spin systems has been an excellent ground for testing ideas of quantum 
information processing in the few-qubit regime. The progression to non- 
trivial numbers of qubits will most likely be accomplished by transitioning 
to engineered single spin systems, and much progress has taken place in this 
direction recently. Promising candidates include quantum dot electron spins 
|1U7] . Nitrogen- vacancy centers in diamond |lU8j . donor electrons/nuclei in 
Si |109] . and Nitrogen atoms trapped in C60 |110] . among others. A com- 
mon theme among many of these approaches is to use the electron spin for 
initialization, fast gate operations and readout, and to use nuclear spins for 
long-time qubit storage or, in the case of quantum dots, as a controllable local 
effective magnetic field jlTI] (the nuclear spin in III-V quantum dots present 
an unfortunate decoherence problem when left uncontrolled |112j ). Single 
qubit control is realized by some form of magnetic resonance in each of these 
approaches (with the exception of the singlet-triplet qubit in quantum dots 
|113j ). whereas two-qubit coupling strategies vary widely. Since these are sin- 
gle quantum systems rather than ensembles, spatial addressing or a mix of 
spatial and frequency addressing becomes possible, an important advantage 
for scalability over frequency-only addressing. These systems are in principle 
more readily scalable than bulk magnetic resonance of molecules, however, 
developing reliable qubits with single quantum systems typically carries a 
host of technological challenges. Though none of the approaches listed above 
has yet moved beyond a small handful of qubits, increasingly higher quality 
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Figure 9: From |106] . (a-d) The average magnetization of the three-spin 
anti-ferromagnetic Ising system versus simulated field h at low temperature, 
showing theoretical predication, numerical simulation of the NMR experi- 
ment including decoherence, and experimental results. Below, entropy 5* 
versus field h, also at low temperature. The entropy peak at /i = indicates 
frustration. 



single and few-qubit systems are being realized at a rapid pace; some ex- 
amples include demonstration of dynamical decoupling |114j and multi-qubit 
control jll5] in NV centers, single-spin readout in Si |116j . and high-fidelity 
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single qubit control and refocussing in quantum dots |117[ 1118] . 

A significant challenge is to achieve fast, high-fidelity readout of a single 
spin (note that projective readout of the spin state is quite different from sim- 
ply detecting the presence of a magnetic moment). So far, the most promising 
methods are optical (e.g. spin-dependent optical transition in NV center) or 
via electron transport (e.g. using the Pauli spin blockade in a double quantum 
dot |113[ 1119] ). Atomic force |120] and nano-magnetometry methods |121j 
are able to detect single spins, but require orders of magnitude improvement 
in sensitivity/ a/ Hz before single-shot spin readout becomes feasible. These 
magnetometry experiments have nonetheless opened up a wide range of pos- 
sible applications in quantum sensing, i.e. exploiting quantum coherence to 
surpass classical limits on measurement sensitivity |121j . and this technology 
may play an eventual role in spatial readout of spin-based quantum proces- 
sors. See |122j for a seminal demonstration of quantum sensing of magnetic 
fields using liquid-state NMR. 

Another common feature of these approaches is the solid-state environ- 
ment surrounding the qubits, which typically leads to shorter decoherence 
times than more isolated systems like ion traps or liquid-state NMR. Al- 
though spin- 1/2 particles are immune to direct coupling with electric fields 
(the dominant noise source in solids), the spin-orbit coupling together with 
phonons provides a pathway for spin relaxation of electrons, which can then 
in turn act as magnetic noise sources for nuclei. Here, dynamical decoupling 
cannot typically improve the situation, because the correlation time of the 
electron spin relaxation process is usually much shorter than the timescale 
of control. Low temperature is thus required in order to reduce the den- 
sity of phonons and suppress relaxation, typically leading to electron spin 
relaxations times in milliseconds for defect centres in dielectric crystals and 
for spins in quantum dots |123j . and up to minutes for electrons at shallow 
donors in high purity Si at 1.2K |124j . The latter work in high purity Si 
demonstrated electron spin coherence times exceeding one second, a ground- 
breaking result for a solid-state system |124j . 

A set of nuclear spins, evolving under nuclear-nuclear dipole couplings, 
can also act a magnetic noise source for an electron coupled to one or more 
of them. This is the case for III-V quantum dots, where an electron spin 
is coupled to ~ 10^ nuclear spins at once due to the contact hyperfine in- 
teraction |125j . The electron spin is dephased on a timescale ~ 10ns due 
to statistical fluctuations of nuclear polarization, but orders of magnitude 
longer dephasing times have been achieved with dynamical decoupling |117 . 
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Similar spin dynamics take place for the electron spin in NV centers which is 
coupled to a small number of proximate natural abundance ^^C nuclei, or for 
donor electrons in Si coupled to ^^Si nuclei; again dynamical decoupling is 
seen to improve coherence times significantly |114t 1126] due to the slow corre- 
lation time of the nuclear bath. In the case of electronic defects in insulating 
or semiconducting crystals, the nearest 'shell' of nuclei may have resolvable 
hyperfine couplings and can therefore be utilized as qubits, whereas more 
distant nuclei have unresolvable couplings and simply generate (dephasing) 
magnetic noise. 

Hybridizing spin and other quantum degrees of freedom, such as pho- 
tons, may solve some of these challenges and is an active area of research. 
The aforementioned NV center is already a type of hybrid system in which 
electron spin can be readout and initialized optically, and furthermore, co- 
herent quantum information stored in the spin state can be converted to 
a photonic 'flying qubit' |127] . The strong coupling cavity quantum elec- 
trodynamics regime has been demonstrated with NV center optical dipole 
transitions coupled to high finesse optical cavities |128l 1129] . opening the 
door to photon-mediated coupling of spatially distant spin qubits. Simi- 
lar efforts are underway to achieve strong coupling between quantum dot 
spins and microwave photon modes in on-chip superconducting resonators, 
for example, using strong spin-orbit coupling to couple the spin qubit to the 
cavity electric field |13Uj . Besides applications in quantum communication, 
this kind of approach could allow for distributing processing tasks optimally 
across different qubit realizations, and for minimizing the number of gates 
needed in a computation by increasing the effective dimensionality of the 
coupled network. 

Chemistry has played a role in the early development of quantum con- 
trol in the context of bulk magnetic resonance of molecular ensembles, and 
one can envision 'bottom-up' architectures emerging for scalable QIP based 
on patterned molecular arrays or monolayers, e.g. those investigated in the 
molecular electronics research area |131l I132j . A viable approach would likely 
use some form of stable radical to provide electron spins for initialization, 
fast manipulation, and readout, and nuclear spins for ancillae and quantum 
memory. The principle challenges in such a system would be addressing and 
readout of spins on individual molecules. Addressing could be achieved using 
suitably strong magnetic field gradients to encode spatial information in the 
frequency domain, and with suitable improvements in sensitivity it might 
be possible to use a nano-diamond NV center scanning probe 'read head' 
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to achieve single-shot spatial readout of spin qubits |121] . The techniques 
for electron-nuclear hyperfine control discussed in section |2.3| would then be 
invaluable tools for implementing a set of high fidelity, universal quantum 
gates in such a system. 
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